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Abstract: 

Evolution equations describe the effect of consecutively integrating out all quantum 
fluctuations with momenta larger than some infrared cutoff scale k. We develop a 
formalism for the introduction of collective degrees of freedom at some intermediate 
scale and derive the corresponding evolution equations. This allows to account for 
the appearance of bound states at some characteristic length scale as, for example, 
the mesons in QCD. The vacuum properties, including condensates of composite 
operators, can be directly infered from the effective action for k — > 0. We compute 
in a simple QCD-motivated model with four-quark interactions the chiral condensate 
ipip and the effective action for pions including A full treatment of QCD along 
these lines seems feasible but still requires substantial work. 
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1 Introduction 



A system with many degrees of freedom or a field theory is often described by 
different relevant excitations at different length scales. A famous example in particle 
physics is the theory of strong interactions, QCD. At short distances (on scales below 
1 fm) the relevant degrees of freedom are quarks and gluons whose interactions are 
well described by perturbative QCD |TJ. The particles which are observed at large 
scales, however, are mesons and hadrons. The interaction of the pseudoscalar mesons 
is well modeled by chiral perturbation theory ||. The corresponding nonlinear a- 
model shares the flavour symmetries of perturbative QCD and it is believed that 
its free phenomenological parameters can ultimately be computed from the QCD 
Lagrangian. In practice, the transition from one set of degrees of freedom (quarks 
and gluons) to another (mesons and hadrons) is a difficult task. Its solution would 
greatly enhance the predictivity of QCD, since the large number of parameters 
characterizing the meson masses and their low energy interactions would be reduced 
to the one free mass scale present in QCD plus the values of the short distance quark 
masses. In this paper we will propose an approach for a description of the "quark- 
meson transition" which is formulated in continuous spacetime. It is sufficiently 
general to apply, in principle, to all problems where the relevant degrees of freedom 
depend on the length scale. 

Our approach is based on the concept of the effective average action T^]. Tk[<p] 
is obtained by integrating out all modes of the quantum field with momenta larger 
than an infrared cutoff scale, q 2 > k 2 . For k — > the infrared cutoff is removed and 
Tfe becomes the generating functional of the 1PI Green functions. In particular the 
properties of the vacuum and of excitations around the vacuum can directly be read 
off from Tfc^o- For nonvanishing k the integration of quantum fluctuations is only 
partial. It does not yet include the effects of modes with q 2 < k 2 . For k — > oo (or k 
equal to some ultraviolet cutoff A) the effective average action equals the classical 
action. Knowledge of the /c-dependence of therefore allows to interpolate from 
the classical action to the effective action, including consecutively more and more 
quantum fluctuations as the scale k is lowered. The scale dependence of the effective 
average action is described by an exact evolution equation^] (t = In k) 

|r tM = iiv((rfM + fl fc )- 1 ^l. (i.i) 



3 In the present form this equation was first obtained || as a renormalization-group improved 
one-loop equation for the average action. It was later proven by simple manipulation of a functional 
integral with an infrared cutoff that eq. (1.1) is an exact nonperturbative equation Q. The 
equivalence of ([bl]) with Polchinski's version of an exact renormalization group equation |5j was 
established by means of a Legendre transform || 0]. (Note that the "ultraviolet cutoff action" 
described by || does not become the effective action as the cutoff goes to zero, but the generating 
functional of one-particle reducible Green functions with sources replaced by specific functionals 
of the fields g, ||, [u^]. A detailed discussion and comparison can be found in pit -) JThe intentions 
of Polchinski as well as subsequent authors || ||, [l^, [b| were mainly aimed towards proofs of 
perturbative renormalizability. Many earlier versions of exact renormalization group or evolution 



equations exist in statistical mechanics |14j. Equation (1.1) is surely equivalent to them. Its 



solutions for k — > are expected to be solutions to the Schwinger-Dyson equation |K 



I 



Here the trace runs typically over momenta and internal indices, and r^ 2 ^ denotes 
the second functional derivative with respect to the fields ip. Within a loopwise 
expansion, = T^ k + T^ k + eq. ( P-.1Q has the iterative solution 

T fc = T (0)k + X - lnDet(rg] fc + R k ) + ... (1.2) 

where the dots denote /c-independent counter terms and higher-order terms. Thus 
equation flLlQ easily reproduces the familiar one-loop result, but the equation is 
actually exact to all loop orders. The additional infrared cutoff term R k in the inverse 
propagator suppresses the propagation of modes with q 2 < k 2 . The appearance of 
T k 2 ^ on the right-hand side turns equation ( |1 . 1|) into a complicated partial differential 
equation for infinitely many variables (typically k and ip(q), the Fourier modes of the 
field ip). Alternatively, equation ( |1.1[) can be viewed as an infinite system of coupled 
nonlinear differential equations for the flow of the infinitely many couplings needed to 
parametrize the most general form of T k . In most cases it is impossible to find exact 
solutions of this equation and its practical use may therefore be questioned. The 
particularly simple "one- loop" form of (|1 . 1|) , the simple interpretation of T k both for 



k — > and k — > oo and the representation of T k in terms of a functional integral with 



constraint JT6[ allow, however, an educated guess on the relevant variables for a given 
theory. This permits the "truncation" of T k with an ansatz containing only a finite 
number of couplings. The existence of useful nonperturbative truncations seems 
to us the main practical advantage of ( |1.1| ) compared to earlier versions of exact 



renormalization group equations [p~4| , ||, whose approximative solutions in general 
involved again usual perturbative, large N or e-expansions. In fact, the formal 
exactness of the evolution equation constitutes by itself not yet any calculational 
progress compared to more standard perturbative expansion techniques. Only the 
existence of nonperturbative approximation schemes for its solutions can provide a 
key for the exploration of nonperturbative physics. 

Despite its close ressemblance to a one-loop equation the evolution equation ( |1.1| ) 
still contains the full nonperturbative content of a theory. We observe that the r.h.s. 
of (|1.1|) is both infrared and ultraviolet finite if an appropriate cutoff R k is chosen, 
as for example 

Rk = q 2 exp(-q 2 /k 2 ) ^ 
1 — exp(—q 2 /k 2 ) 

lim R k = k 2 (1.4) 

No need for an additional specification of an ultraviolet regularization arises in this 
case. This is replaced by the necessary specification of T k for some sufficiently 
large scale k = A. (Here A may be associated with a physical ultraviolet cutoff 
in the sense that all quantum fluctuations with q 2 > A 2 are already included in 
Ta-) A solution of the exact evolution equation for k — > with initial condition 
specified by the form of T\ yields all Green functions of the theory. The possibility 
to extract useful nonperturbative information even from approximative solutions of 
eq. (|1.1| ) has already been demonstrated by the successful use of nonperturbative 
"truncations": For a SO(N) symmetric scalar field theory the phase structure in 
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arbitrary dimension d has been correctly reproduced |J, including the Kosterlitz- 
Thouless phase transition for d = 2, N = 2. Critical exponents of the three- 



dimensional theory have been computed with a few percent accuracy |L8[ . The high 
temperature phase transition in the four-dimensional theory was established to be 
of the second order [Tj| in contrast to indications of earlier perturbative results. 
(For more related nonperturbative work on scalar field theories see refs. [2C|| .) Of 
direct relevance to the present work is the description of bound states @ through 
properties of T^: A nonperturbative study of the momentum dependence of the four- 
point function shows the appearance of a pole in the s-channel from which mass and 
wave function of the bound state can be extracted. These first encouraging results 
for the description of nonperturbative physics raise the challenge: Can one tackle 
with these methods the most important nonperturbative problem in field theory - 
the behaviour of QCD for low momenta? 

One problem for such a project is immediately apparent: For large k (k ^> 
1 GeV) the theory is well described in terms of quarks and gluons and Tj. should 
take a relatively simple form in terms of these fields — typically dominated by 
gauge-invariant kinetic terms and mass terms for the quarks. Due to the growing 
gauge coupling one expects that the form of Tk gets more and more complicated 
as k comes close to typical scales where strong interaction phenomena manifest 
themselves, say at k around 1 GeV. At these scales effective four-quark couplings or 
terms ~ (F^F^) 2 do no longer remain small corrections. For k even much smaller 
than 1 GeV the dominant excitations are mesons and hadrons. A description of 
the meson and hadron physics in terms of the original quark and gluon degrees 
of freedom must be extremely complicated: Bound states result in a complicated 
momentum dependence of the n-point functions and condensates indicate that an 
expansion of around the configuration of vanishing fields presumably becomes 
meaningless. One clearly needs a transition to express Tk in terms of new degrees of 
freedom more adapted to the physics. In this paper we want to provide the necessary 
formalism to include in collective fields (as appropriate for the mesons and ipip 
condensates). We will again derive an exact evolution equation equivalent to (jl.lQ , 
but now including the collective fields. Even though on the formal level there is no 
need to use these new "collective field equations" the necessity of an appropriate 
nonperturbative truncation scheme makes such an approach almost compulsary. The 
present work builds on earlier descriptions of collective fields [[Kj in the framework 
of Polchinski's exact renormalization group equations 

To be more specific, we expand T k [ip] in powers of the "fundamental" field ip, 
with fc-dependent one-particle irreducible n point functions 1^ as coefficients. (In 
the QCD example ip stands for quark and gluon fields.) The n-point functions 
rj^ will possibly have a complicated dependence on the external momenta. Exact 
flow equations for each of them can be derived by taking appropriate functional 
derivatives of (|ITT|). Let us, for simplicity, consider two-particle bound states or 
condensates in the following: Informations on possible two-particle bound states 
are contained in the four-point function T^ 4 ^ and are usually extracted by means 



of the Bethe-Salpeter equation plfl . In our approach we follow the /c-dependent 
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four-point function T k from large towards smaller values of k. If the theory under 
consideration possesses a two-particle bound state in the channel (f\(f2 — >■ ^3^4, ll 
will develop a pole in the momentum variable s = (pi +P2) 2 for small enough 
In this case T k ^ assumes the form 

rl 4) (pi, ...Pa) ~ -^(Pi,P2)G , (s)^(p 3 ,p 4 ) + rj 4) (pi, ...,p 4 ) (1.5) 

where g(pi,P2) corresponds to the (amputated) Bethe-Salpeter wave function and 
G(s) to the bound state propagator approaching a pole-like dependence on s for 
small enough k. At the pole rj 4) (pi) remains a bounded function of the momenta 
Pi, ...,p&. In ref. 0, where the suitably truncated flow equation for 1^ was inte- 
grated numerically for different models, the emergence of the structure described in 
( p..5|) was indeed observed. It was demonstrated that both the wave function g and 
the propagator G(s) near the pole can be reliably extracted from rjj. . 
We associate the bound state with a composite operator of the type 

O(q) = J ^^g{pi,P2)^\-pi)^P2){2^)H\q-p 1 -p 2 ) (1.6) 

A collective field a is then introduced for the composite operator O. (In our QCD 
model <p stands for the quarks and a denotes the mesons of the linear a- model.) 
Since the collective field a associated with the bound state may also develop a 
nonvanishing vacuum expectation value, it would be desirable to study the complete 
effective potential of the field a and not just its part quadratic in a. (The latter 
could be obtained from G(s) at s — 0.) This is one of the main motivations for 
the formalism to be developed in the next sections. This formalism deals with 
the effective average action T^Up, a] involving one (or more) collective fields a. It 
requires, at a certain scale k v , to replace Tfc^f^] by (cf. eq. ( |5.11|) ) 

T fc >,a] = r fc >] + l ~0%]GOUp] - JO[<p] + l^G-'a (1.7) 

The complete effective action for both fundamental and collective fields TUp, a] is 
then obtained after evolving Tk[(f, &} from k = k v down to k = with the help of a 
modified flow equation discussed in sects. 3, 4. 

If the four-point function rjj. has, at a scale k ~ k v , developed a pole-like 
structure as in (|1.5| ), a corresponding term will appear in the part quartic in tp in 

Up] on the r.h.s. of (|1.7|) . The operator OUp] can then be chosen appropriately 
such that the pole-like structure in T k ^ is cancelled by the second term on the r.h.s. 
of ( |1.7| ). With this choice the operator OUp] contains the two-particle wave function 
of g(pi,P2) of eq. ( |1.5|) . After the cancellation of the pole-like structure in 1^ 
the evolution of r^f^, a] down to k = can proceed with new "initial condition" 
~ 1"! ■ The remaining part of the four point function stays then bounded 

4 Even though we work in Euclidean space where the evolution equation involves with all 
Pi>0 and s > 0, the final expressions for the Green functions can be analytically continued to 
Minkowski space and s < 0. 



4 



(unless further bound states appear). In this way it becomes apparent how the 
introduction of the collective field a can remove a dominant strong interaction with 
a complicated momentum dependence from the Green functions of the fundamental 
fields. The price to pay is an effective average action with more fields and therefore 
more possible couplings. Even though at the scale k v the effective average action 
starts only with terms linear and quadratic in the collective field a as in eq. ( |1.7| ), this 
property will not be preserved for smaller values of k. The flow equation typically 
generates interaction terms for a as well. Despite the apparent complication through 
the presence of additional fields, the approximate (truncated) flow equations may 
become much simpler. Since we have now explicitly included the relevant degrees of 
freedom, the important physics may be described in terms of fewer couplings than 
needed for a reliable description of in terms of the fundamental fields alone! 

This paper is organized as follows: In sect. 2 we develop the formalism of the 
effective average action including collective fields, and we derive exact flow equations. 
In sect. 3 we present the flow equations in a version where "fundamental" and 
collective fields are treated on equal footing. Various strategies how to add or 
remove effective degrees of freedom in dependence on the characteristic length scale 
of the problem are discussed in sect. 4. For the case of QCD this gives the formalism 
how to replace quark degrees of freedom by meson degrees of freedom. In sects. 5 
and 6 we finally apply our formalism to a QCD motivated model with four-quark 
interactions. We integrate out the quarks and obtain an effective potential for the 
linear cr-model. This allows a determination of f n and the chiral condensate ipip. 
Conclusions and prospects for quantitative improvements of our results are contained 
in sect. 7. 



2 Exact evolution equations with composite op- 
erators and collective fields 

Our starting point is the scale-dependent generating functional for the connected 
Green functions Q 

W kfi [J,K\ = In J Vxexp-Sklx] (2.1) 

S k [x] = S[ X ] + A k S[ X ] - Jh - K^GO[ X ] - \k^GK (2.2) 

Here x a are the (bosonic) quantum fields of the theoryf]with action S[x] and J* the 
associated sources with 

J ] X = J* a X a (2-3) 

The indices a include momenta q 11 , possible Lorentz indices for vector and tensor 
fields as well as internal indices. (For example, the action for free scalar fields 
in d Euclidean dimensions reads S = \x^<fx — §/ p^Xa(?)(?M5 M )x°(9)- m this 

5 The second index of Wf-o refers to an optional infrared cutoff for collective fields to be 
discussed in sect. 3. 

6 The generalization for fermionic fields is straightforward p2[ and will be given in Appendix C. 
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notation functional derivatives obey Sx a (q)/Sx b (q') — ^(9> where S(q,q') = 
(2ir) d 5 d (q - q 1 ) such that Tr q S(q, q') = 1.) 

Let us neglect first the last two terms in Sk which involve K. Then the only mod- 
ification of the standard definition of the generating function W[J] is the addition 
of an infrared cutoff which is quadratic in \\ 

&kS[ X ] = \x ] RkX (2.4) 

The matrix (Rk)p will be chosen such that it suppresses the contribution of fluctu- 
ations with momenta q 2 <C k 2 in the functional integral ( |2.1|) . We require that the 
infrared cutoff Aj,S preserves all symmetries of the theory. A typical choice 7 for 
for the case of scalar fields with internal indices a, h is 

(Rk)\(q,q') = R t (qK5(q,q') (2.5) 

W) = f^lfi (^) 
1 — exp(— q 2 / k 2 ) 

For q 2 <C k 2 one finds Rk ~ k 2 and the infrared cutoff scale k acts like an additional 
mass. On the other side Rk is exponentially small for q 2 ^> k 2 and does therefore 
not affect the functional integration of modes with q 2 ^> k 2 . If desired, however, 
Rk(q 2 ) could also be chosen such that it incorporates an UV cutoff A. Then only 
fluctuations with q 2 < A 2 are included in the functional integral (see sect. 5). A 
very intuitive picture arises if we combine the classical kinetic term from S[x] with 
AkS which leads to a lowest order propagator G k °\q) of the formQ 

4%) = tf + ^-'= 1 ~ exp( r g2A2) (") 

q z 

We see that the propagation of modes with q 2 3> k 2 remains essentially unchanged 
whereas modes with q 2 k 2 cease to propagate. 

We observe that for k — > and K — the functional Wk t o approaches the 
standard definition of W[J) 

lim A k S = 

fc->0 

\imW k>0 [J,0}=W[J] (2.8) 

Performing a Legendre transformation of Wk,o for K = and subtracting the infrared 
cutoff term yields the effective average action T^qM which becomes the generating 
functional for the 1PI Green functions for k — ► M 

r fc) oM = -w k ,o[j, o] + jV - 



7 For precision calculations it is convenient to incorporate in Rk(q) appropriate wave function 
renormalization constants and possibly also mass terms. 

8 If Rk{q 2 ) involves an UV cutoff A, it can be chosen such that Gjj^\ has the form — 
({exp-q 2 /A 2 )-exp(-q 2 /k 2 ))/q 2 . 
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= Ml (2.9) 

5 J* 

The dependence of I\ on the scale k is described by the exact evolution equation 
( [LID mentioned in the introduction. 

We now want to consider composite operators % [x] for which we have introduced 
appropriate sources K* in (|2.1|) , with 

K^GO[x} = K*&^[ X ] (2.10) 

Here G typically plays the role of a propagator with an appropriate dimensional 
factor such that the sources have a standard normalization. In momentum space G 
is a function of q 2 . For simplicity of notation it is often convenient to use rescaled 
operators 

0[x] = GO[ X ) (2.11) 
such that the coupling to the sources takes the standard form K^O. We also have 



used the freedom |23], ^TJ to add in (|2T2| ) a field-independent term quadratic in K 
with a matrix & y The role of this term will become more apparent later. The 
generating functional for the connected Green functions for the fields x as we U as 
for the composite operator obtains again in the limit k — > 

\imW kfi [J,K} = W[J,K] (2.12) 

For example, the expectation value of the operator % reads 

SW 

< O l >= ^- (2.13) 



UJX l \J=QK=Q 



and the connected two-point function isf] 



< &0) >-<&>< 6* >= - & j (2.14) 



|j=o,e-=o 



Knowledge of the functional W[J, K] therefore allows the computation of various 
"condensates" and correlations between "condensates" . 

More generally, we will define for an arbitrary scale k and arbitrary source J, K 
classical fields <p a and a\ 

SW k ,o 



5K* 



a\ (2.15) 



9 We have chosen a normalization appropriate for real composite operators where O* and Oi 
are not independent. Our conventions also cover complex fields and operators if index summations 
over internal indices include negative values 
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where the composite or collective field a % represents the composite operator 0{. We 
combine the fundamental and composite fields and their associated sources in single 
vectors 

3m = i^ai ) 



ij m (2.16) 



Si* 

The generalized two-point function reads then 

G m n = (2.17) 

We now perform a Legendre transform with respect to both "fundamental" and 
"composite" sources J and K 

f M N + ^,o[j]-jV = (2.18) 

8Y k n 



Sijj r 



fra (2-19) 



G ™tPtt = 5 v (2- 2 °) 

In particular, the values of if and a extremizing I^o determine the expectation 
values of x an d 0[x] an d the second functional derivative of I^o at the minimum 
gives the exact inverse propagator. 

The exact evolution equation or flow equation describing the dependence of i\ 
on the infrared cutoff scale k is now easily obtained. We notice that the only 
dependence of W kj o on k arises through A k S which is quadratic in x- The derivative 
of Wk,o with respect to t = In A; can therefore be expressed in terms of the Green 
function G 

|r M| , = —w^ = | < A k s > 



id_ 

2di 



(R k ) a p < X* a X P > 



K^)/^ +</><x>) (2 - 2i) 



Subtracting from f ko the infrared cutoff 



\^Rk<P (2-22) 



and using the identity ( ^20|) gives the final form of the exact evolution equation for 



the effective average action T k including composite fields 
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It expresses the scale dependence of I\o in terms of its second functional derivative 



( 



1-1(2) \ m ^ 2 Tfc 

w n = s^sr (2 ' 24) 



We have extended in (|2.23| ) the matrix Rk to act on vectors ■0, with (Rk) l j = 
0, {Rk) 1 a = 0, (Rk) a i = 0. Only the (a, (3) components of Rk (the ones corresponding 
to the fundamental fields p a ) depend on k in this formulation. Even though only the 



a,P components of (T kQ + Rk) 1 contribute therefore in ( |2.23|) , the (i,j) and (i, a) 

components of T k 2 Q are relevant since they enter in forming the inverse of T^q + Rk- 
For k — > oo i\ [y3,cr] approaches 

rW^a] = S[<p] + ±Oi[<p]GO[<p]-JO[<p] + i^G-V. (2.25) 

In a theory with an ultraviolet cutoff A (cf. footnote 7), the average action becomes 
equal to the r.h.s. of eq. fl2.25| ) for k = A. Solving the evolution equation fl2.23|) 
with this initial condition interpolates from the short distance physics described by 
T A to the physics at longer distances described by T k , k < A. In this process the 
quantum fluctuations with momenta A 2 > q 2 > k 2 are integrated out. Tk for k — > 
therefore amounts to the complete effective action of the quantum field theory. 

In summary, we have presented here an extended form of the exact evolution 
equation (|1.1[) . It uses additional composite fields a. JA solution of Tk[<fi, cr) for k — > 
encodes all information of the original procedure, namely all 1PI Green functions 
for the fundamental fields x- m addition, it contains the complete information on 
expectation values of composite operators 0[x] and their correlations. At first sight 

(2) 

the use of an enlarged matrix T k ' may seem to be an additional complication which 
is the price for gaining additional information about Green functions for composite 
operators. It may happen, however, that physically important pieces in Tk have 
a rather complicated form when expressed in terms of the "fundamental fields" 
(p, but become simple once they are written in terms of composite operators. In 
this case the use of more fields (ip, a) may result in a considerable simplification 
of calculations. This applies in particular to models with condensates of collective 
fields and/or bound states. For the theory of strong interactions we may associate tp 
with the gluon and quark degrees of freedom and a with operators like 'ijjijj, F IMV F* lv 
or mesons, hadrons, and glueballs. The low energy part of the effective action is 
expected to become much simpler when expressed in terms of a rather than of p. 



3 The two-field formalism 

In a variety of physical situations there is no clear distinction between fundamental 
and composite fields. This applies, for example, if one deals with processes with 
typical momentum transfers much smaller than the mass scale charcteristic for the 
structure of the bound state. The structure cannot be resolved in this case, and 
the bound state appears as an ordinary particle. An example are versions of the 
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standard model where the Higgs scalar is a bound state. In situations of this type 
one would like to treat fundamental fields and composite fields on a completely 
equal footing (at least at low energies) such that standard field-theoretical methods 
as perturbation theory can be applied. For this purpose it is instructive to write the 
functional W kt o[J, K] in terms of a functional integration over both fundamental and 
collective degrees of freedom. We use the identity (choosing G such that & = G) 

l = Njv P ex P {-I(pt - tfG - 0^ X ]G)G-\p - GK - GO[ X })] (3.1) 



with N an irrelevant constant to be omitted in the following. Insertion in (2.1) 
yields 

W k , [J,K] = In/ D X Dpexp-S k [ X ,p} (3.2) 

with 

S k [x,p} = S[ X ]-Jh + ^kS[x] 

+ l -0^ X }GO\x] - p ] 0\ X ] - K^p + l -p^G- l p. (3.3) 

We see that Wj^of./, K] has now the same form of a theory with additional fields and 
interactions. In particular, S k is linear in the source K which motivates the term 
~ K^GK in the original formulation ( |2.2|) . Compared to the formulation in terms of 
only fundamental fields x we have introduced new fields p with inverse propagator 
G^ 1 and interactions with the fundamental field ~ p^O[x]- For the example of an 
operator O quadratic in x one obtains a cubic interaction as depicted in fig. 1. One 
also has an additional term ~ O^GO involving only the fields x- It h &s the form of 
an exchange of p in the tree approximation (fig. 2) and can be used to cancel terms 
in S[x] with a pole structure as discussed in the introduction and in more detail in 
sections 4 and 5. 

Inside the one-loop diagram corresponding to the r.h.s. of the flow equation 
( |2.23 ) also collective fields appear. Their presence is induced by the inversion of 



(2) 

(r k ' Q + R k ) in the enlarged space with indices (i, at), (cf. the discussion following eq. 
( j2.24| )). The effective propagator G of the collective fields p does not yet include 
an infrared regulator; thus it is possible that the r.h.s. of the flow equation fl2.23|) 
contains, in a given truncation, an infrared divergence. This situation can especially 
arise in the case of dynamical symmetry breaking, where some collective degrees of 
freedom become massless Goldstone bosons]^ In the formulation Q3.3p it becomes 
easy to introduce an additional infrared cutoff (with scale k) for the collective field 
p. We can generalize W k ,o to W k ~ k by adding to S k ( |3.3|) a term 

A k S[p] = ^R kP . (3.4) 



10 In practice one might sometimes consider a truncation of T k where vertices generating "in- 
ternal" collective fields are discarded, see sect. 5. Here, however, we consider the general case. 



10 



The matrix (R k Yj should have similar properties as (Rk) a p. In particular it should 
vanish for k — > 0, and all eigenvalues should diverge for k — > oo (or k — *■ A). 
Correspondingly, T k k is defined by the analogue of (p. 181) and T k k reads 



We note 



1 1 

T k,k = T k,k ~ I}? R W ~ ^R k o. (3.5) 



with the expectation value evaluated now with the action Sfc[x, p] + A k S[p\. The 
evolution equation for the ^-dependence of Y k k at fixed k has the same form as ( p.23|) 

with (Rk) l j replaced by (R k Y which is now different from zero but independent of 
k. (Only the denominator on the r.h.s. of (p.23|) changes due to the new definition 
Q3.5|). As in the previous section, (Rk) a /3 and (R k Yj are combined to an enlarged 
matrix R kk .) The dependence of r fc k on t = In k can also be derived in complete 
analogy to ( ggg ) 

4r a = iTv(^i(rS + Bi ,)-'l (3.7) 



Here £ acts on vectors (y?, a) as defined above and only -§%{R kk Yj = ~§t( R k Yj * s 
different from zero in this case. In particular, we may identify k and k and define 

T k [ip,a] = F kjk [ip,a] (3.8) 

The evolution equation for the /c-dependence of T k is again given by Q2.23| ) since 



d d d 

~dlF k = dt^ k M~k=k + Q^k,k\k=k (3-9) 

Now both (R k ) a f3 and (R k Yj depend on k. The evolution equation for the two sorts 
of fields if and a has exactly the same form as the original equation ( |1 . 1| ) . 

In the limit k — > A, k — > A (A may be infinity) the functional integrals defining 
W k = W k ,k and Y k are easily solved. The term quadratic in the fields ~ R k diverges, 
the classical approximation becomes exact and one obtains the same result as in 
(PI): 

V A [<p, a) = S[<p] + hj%\GO[ V ] - JO[<p] + X -^G- X a (3.10) 

The equivalence of the initial condition ( |3.10| ) with the original formulation in terms 
of fundamental fields only is readily established by solving the field equations for a 
and inserting in T A 

da \a 

a = GO[(p] = oy\ 

r A [<p,a ]=S[tp] (3.11) 
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In the opposite limit k — > 0, k — > the infrared cutoff term ~ Rk vanishes and 
To becomes T^o of sect. 2, the generating functional for the 1PI Green functions for 
the fields tp and a. Inserting the solution of the field equation for a which minimizes 

^ =0 (3.12) 

da \ ao 

we recover the generating functional for the 1PI Green functions for the fundamental 
field x (since K = 0) 

r[<p] = r [<p,(T [tp]] (3.13) 

The full generating functional Fo[y?, cr] contains all information on the vacuum ex- 
pectation of the operator 0[x] as well as all n-point functions of this operator. The 
procedure for extracting these quantities is straightforward and given explicitly in 
appendix A. 



4 Scale-dependent degrees of freedom 

The relevant degrees of freedom often depend on the length scale. For the example 
of QCD one would like to describe the short-distance physics in terms of gluons 
and quarks and the long-distance physics in terms of mesons and hadrons. In sect. 
3 we have introduced the effective average action r fc ~ k depending on two scales k 

and k, which will allow us to describe a smooth transition between different sorts 
of relevant degrees of freedom. Consider first the case that below some momentum 
scale k a we want to describe the physics uniquely in terms of the collective fields a. 
For example, one may want to describe QCD below k a ~ 500 MeV by a linear or 
nonlinear a-model for the pseudoscalar mesons. In our approach this can be done 
very naturally by first fixing the infrared regulator for the collective fields k = k a 
and solving the evolution equation for the infrared regulator for the fundamental 
fields k — > 0. The resulting effective action for composite fields T^fcr] is given by 

riM = r^M*U] (4.i) 

where 

^[<PoW°] = (4-2) 

for k = k a . It can be used as a starting point for integrating out the "mesons" 
with momenta q 2 < k 2 . The low momentum fluctuations of the mesons are not yet 
accounted for in T^[a] due to the existence of an effective infrared cutoff ( |3.4| ). They 
will be included by evolving T^a] to T [a] using the evolution equation ( |3.7|) . We 
observe that R is now nonvanishing only for the composite part {R a p = since k = 0, 

R l j = (Rj^Yj) and that T^ 2 - becomes block-diagonal in composite and fundamental 
fields 

S 2 T - 

-^-^ = (4.3) 

6a6i P \<fo,<r 
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provided <po[o] ( |5.2| ) is independent of a. (If the fundamental fields ip assume no 
vacuum expectation values, the solutions of eq. (4.2) will simply be <po[cr] = 0.) 
Then the evolution equation ( |3.7]) only involves r~ [a] on the r.h.s. and has exactly 
the same form as ( |1.1|) , now expressed in terms of composite fields. This statement 
generalizes to an arbitrary dependence of ipo on a as briefly explained in appendix 
B. We conclude that the solution of T k ^ for k — > at fixed k can be used for a 
complete replacement of the fundamental variables ip by the composite variables a. 
This procedure amounts to integrating out the fluctuations y? in a background 
The same procedure can actually also be applied to the problem how to integrate 
out only one sort of fundamental fields (without reference to collective degrees of 
freedom). One may introduce a variable infrared cutoff with scale k only for the 
fields to be integrated out, while the remaining fields have a cutoff with scale A. 
Lowering k to zero at fixed A then only includes quantum fluctuations of the fields 
to be integrated out, and ( |4.1|) , ( |4.2j ) with k = A give the effective action for the 
remaining fields a. 

Under certain circumstances it may be convenient to consider both fundamental 
and composite degrees of freedom on an equal footing. This is done most easily 
by a study of the evolution equation ( |3.9| ) where k = k. An example are top con- 
densate models at scales much below the composite scale [p4], which have already 
been treated by methods closely related to the present one P5fl . Treating the scalar 
bound state (the composite Higgs doublet) similar to the fundamental fermions will 
reproduce most easily the results of the perturbative standard model since for k = k 
the fermion and scalar loops described by the formal expression (3.9) are treated in 
the same way. This is not guaranteed in the ("asymmetric") version (2.23) where 
k = 0. Since the standard contribution of the scalar loops is missing, a relatively 
complicated truncation of I\o may be needed in order to reproduce even the stan- 
dard one-loop result for the running of weak couplings. Another typical application 
of the running of is the transition region between the quark description and 
meson description in QCD at scales between 500 MeV and 2 GeV. 

Bound states are, in general, characterized by a typical "composite" scale. At 
length scales much shorter than the inverse composite scale, they play no partic- 
ular role. For large enough k not only the introduction of a second scale k seems 
cumbersome, but the whole formalism with composite degrees of freedom seems 
sometimes not to be very well adapted. One does not want to describe asymptotic 
short-distance QCD by dealing explicitly with mesons! This problem can be coped 
with by "switching on" the collective degrees of freedom only for the evolution at 
k smaller than some scale k v . This can be achieved by an appropriate choice of 
the composite operator 0[x]- Actually, already the presence of a wave function 
g{PiiP2) m 0[<p\ (see eq.(1.6) ), which decays for large momenta p 2 , will suppress 
the contributions of the new terms to the flow of Tk[tp, a] for large k 2 . This effect 
is much enhanced by the introduction of exponentially decaying functions of pf into 
the operator 0[ip\. Let us represent O as a functional of the Fourier modes x(<?) in 

n The actual calculations of the present paper in sects. 5,6 will not involve internal propagators 
of the CT-field, thus the introduction of a scale k a is not necessary. 
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the form 

0[ X ] = 0[f kv (q) X (q)} (4.4) 
/^(g) = exp(-g 2 /2^) (4.5) 

Typically O may be some polynomial of fk v (q)x('l) involving appropriate wave func- 
tions — see the next section. More general forms of O are also allowed, provided it 
does not contain an exponentially growing momentum dependence. The appearance 
of the exponential suppression factor fu in the coupling between p and x (H-D w iU 
then effectively "switch off" all effects of the composite fields for q 2 ^> k 2 ^. Indeed, 
for k > Ak v with A sufficiently large (say A = 10) all contributions from composite 
fields in the evolution equation ( |2.23| ) are suppressed. This allows us to compute 
rAfc„[^] by integrating the evolution equations for Ak v < k < A without using the 
formalism with composite fields (solving ( |1 . 1| ) ) . On the other hand, the modified 
operator 0[ip] still describes appropriately the "bound state" for momenta p 2 < k 2 ^ 
of the fundamental fields. At k = Ak^ one may switch to the two-field formalism. 
The initial terms on the r.h.s. of eqs. ( |2.25|) or ( p.lOj ), which formally have to be 
introduced for k — > oo (or k — A), did not affect the evolution of T k for k > Ak v , 
and can thus be introduced at k = Ak v only: 

r AK ,M = ?AkM + \<$[<p\GO[<p\ - JO[<p] + Vg-V (4.6) 

At this stage one has several options for the infrared cutoff k for the collective fields 
a, as described before: One can put it equal to zero right away, or one can put 
it equal to Ak^ and treat it as a variable independent of k subsequently (using 
eq. (2.23)), or one can identify it with k both at the starting point k = Ak v and 
concerning the subsequent evolution, using eq. (3.9). Note that the r.h.s. of eq. 
( |16|) does not depend on the choice of k; the dependence of F Ak ^ on k manifests 
itselfe only in the form of the evolution equations. The collective degrees of freedom 
will now affect the evolution of T kk -[<p, a] for k < Ak v due to the disappearance of 
the exponential suppression factor ( [4.5| ) in the a — ip coupling. 



The concrete choice of the operators 0[(p] and the form of the propagator G 
for the collective fields a depends on the problem under consideration. For the 
case of propagating bound states it has been sketched in the introduction. Let 
us concentrate in the following on a "two-particle bound state" showing up in the 
four-point function r^ 4 . Our starting point is the effective action Tk[<p] without 
collective fields a. After the integration of the corresponding flow equations ( |1 . 1|) 
down to a certain "bound state" scale it is assumed that a pole-like structure as 
in (|1.5|) has emerged as part of the four-point function rl . It is now desirable to 



k 

choose the operator 0[(p] and the propagator G such that this pole-like structure 
within gets cancelled; then the collective field a associated with the operator 
0[ip] corresponds to the bound state degree of freedom generated by the dynamics of 
the theory. (Formally no error is made if 0[(p] is not chosen appropriately, but then 
informations on the bound state are still contained partially in the remaining part 
of the four-point function of the fundamental fields.) An obvious possible choice is 
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to identify the ansatz for G(q 2 ) in eqs. (|2.2j), (3.1) and (4.6) with the dynamically 



generated G(s) in eq. ( |1.5|) and to take the Fourier components of 0[tp] as 

O(q) = J ^^/^(pl)/^(P2)^(pi,P2)^(pi)^(p 2 )(27r) 4 5 4 (p 1 + p 2 - q) (4.7) 

Here the scale k v should be somewhat above the scale where a pole-like structure 
in emerges. This ensures that fu does not switch off the collective fields in 
the momentum range where they are needed. On the other hand the presence of 
fu justifies the introduction of the additional terms in eq. ( |4.6| ) only at k — Ak v 
and not already at k —>■ oo (or k = A) as in eqs. (2.25), (3.10). In practice it 
might also be useful to try modified forms of the operator 0[tp] and G introduced at 
k = Ak v ; the optimal choice is the one which completely prevents the appearance 
of a pole-like structure in rjj. for k <C k v . 

The methods discussed in this section allow for describing smoothly the transition 
from perturbative QCD to the nonlinear cr-model. Variations of the transition scale 
Ak v and the different options for the infrared cutoff k for the collective fields would 
allow for a variety of checks. Obviously, much work needs to be done before such a 
program can be implemented. We nevertheless will demonstrate the viability of our 
method by the investigation of a simplified model in sect. 5. We finally mention 
that the restriction of the formalism developed in the last sections to scalar fields is 
by no means necessary. The generalization for the inclusion of (chiral) fermions is 
presented in appendix C. Inclusion of gauge fields requires more thought but seems 
not to encounter insurmountable difficulties M, O, p6 . 



5 A QCD-motivated model 

So far we have developed the general formalism for the description of bound states 
and collective fields in the context of the average action in a very general but nec- 
essarily also rather abstract way. In this section we want to give an example which 
shows how these ideas work in practice. Our aim is a demonstration of feasibility 
of our program rather than an attempt to compute with precision. We nevertheless 
take a model which is inspired by QCD. It contains sufficiently many details such 
that it can account for mesons and chiral condensates and can later be extended to 
full QCD (see the discussion in sect. 7). We work with massless quarks even though 
the introduction of quark mass terms poses no particular technical problem. 

We start with a theory with only fermion fields ip. They correspond to massless 
quarks and play the role of the fundamental fields <p of the last sections. We have 
in mind a version of QCD, where the gluons have already been integrated out. As 
a result we expect that, among others, a (generally nonlocal) four-quark interaction 
T^\pi) has been generated at a scale A ~ 1.5 GeV. Concerning the form of Ia (pi) 
we assume that for pf < A 2 it depends just on the Mandelstam variable t such that 
it corresponds to the sum of a linear and a confining potential in the nonrelativistic 
limit. We consider an expansion of the effective action up to fourth order in the 
quark fields, and we neglect the running of the effective quark propagator. (Due 



15 



to chiral symmetry the quarks will remain massless anyhow; we only neglect the 
running of the quark wave function normalization.) Thus we just have to integrate 
the evolution equation for 1^ with the above 1^ as boundary condition at k — A. 
Now we observe indeed that 1^ at some small scale ~ develops a pole-like 
structure as in (1.5). It allows us to read off the wave function g{jpx,P2) and the 
collective field propagator G(s). 

Next we neglect the difference between the bound-state scale, k v and Ak 9 , and 
construct T[ip,ijj,a] at the scale k v according to the rule eq. (4.6). (Here we put 
fk = 1 for simplicity.) By this procedure the model is transformed into a linear 
cr-model with momentum dependent Yukawa couplings between mesons and quarks. 
The part of involving the scalars starts at this scale only with a quadratic term 
and Yukawa couplings as given by the r.h.s. of eq. (4.6). At scales k below k v we 
focus our attention on the effective potential for the collective field a. It will be gen- 
erated after integrating out the quark fields if), if> with momenta q with < q 2 < k 2 ^, 
due to the (momentum dependent) Yukawa coupling. The dominant effect is the 
"quadratic running" Jig, g2| of the scalar mass term. The fermion fluctuations 
induce a negative mass term at the origin (a = 0) and trigger therefore the spon- 
taneous breaking of chiral symmetry. We neglect contributions to the effective po- 
tential due to internal cx-lines. This allows to compute the effective potential from 
a simple "one-loop formula" (eq. (|6.19|) below) in terms of the parameters present 
in T kv [ip,ifj, a]. 

To be more precise, we describe the fundamental quark fields by four component 
Dirac spinors if)^,if) l a wnere t ne flavour index a runs from l...Nf, i = l...N c is the 
colour index and we omit the spinor indices. The effective action contains terms 
quadratic and quartic in the quark fields, and is taken to be invariant under SU(N C ) 
and chiral U(Nf)y <g> U(Nf) A symmetries: 

5 v if) a = i6 Vb ip\ 5 A r = ils6 Ab ilj h 

5v4>a = -i4> b b Va , S A ip a = iip b 9 b Aa j 5 . (5.1) 

The quadratic part describes the inverse quark propagator, and a possible k de- 
pendent wave function normalization is put equal to 1 (we contract over the not 
explicitly written spinor indices.): 

r 2lfc = / -^-M^mi) (5.2) 

Even in the presence of the vector and axial vector symmetries (5.1) a number of 
different four-quark interactions can be written down. We restrict ourselves to one 
particular spin, flavour and colour structure: 

r 4;fc = 7j/ll ((^j) (^^(P 1 +P2 ~ P3 - P4)Afc(pi,P2,P3,P4) (5.3) 

This structure colour singlets with nontrivial flavour in the s-channel (they will later 
correspond to the mesons) and flavour singlets transforming as colour singlets and 
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octets in the t-channel (as mediated, for example, by a one-gluon exchange). As 
usual the kinematic variables s and t are defined here by 

s = (p x +p 2 ) 2 = (p 3 +Pif 

t = (Pi -Pa) 2 = (pa -P4) 2 - (5.4) 

The contractions of spinor indices are indicated by brackets, and the invariance 
under (5.1) is readily checked. 

The sum of r 2i fc and T^ k has now to be inserted into the fermionic version of the 
flow equation: 

d t T k = -Tr^(T^ + R kF ry (5.5) 

1 - exp(-q 2 /A 2 ) + exp(-q 2 /k 2 ) 
RkF[q) ~ i exp{-q 2 /k 2 )-exp{- q ^B) (5 ' 6) 
(2) 

such that the sum of the terms T k + R kF gives 



We use hereQ 



exp(—q 2 /A 2 ) — exp(—q 2 /k 2 ) 



CW). (5-7) 



We are interested in the terms ~ ('0V') 2 of eq. (5.5), which contribute to the running 
of the coefficient \ k (Pi, P2-> Ps-> Pi) °f After a corresponding expansion of the 

r.h.s. of eq. (5.5) the equation describing the running of X k assumes the schematic 
form shown in fig. 3. Here the lower inner line denotes the regularized fermionic 
propagator as derived from the inverse of rf } + R k , 

G 2 k ( Q ) = ^p(~<l 2 /A 2 )- exp(- q 2 /k 2 )^ 

i 

and the upper crossed line its derivative with respect to t — In k. (In fig. 3 we have 
not shown a similar diagram with the cross on the lower line.) Actually, inserting 
the ansatz (5.3) for T^ k into the r.h.s. of eq. ( |5.5|) , also different spin, flavour, 
and colour structures than the one of eq. Q5.3p are generated on the l.h.s. of eq. 
( |5.5| ). They are put to zero in our truncation. Since only the retained contribution 
is proportional to a combinatorial colour factor iV c our truncated evolution equation 
becomes exact in the leading order in a 1/N C expansion. After taking the various 
combinatorial factors into account and performing the trace over spinor indices, the 
flow equation for becomes 

d t Xk(Pl,P2,P3,P4) = 

— j?rf pSjrA fc (pi,P3, Q, ~q + Pi+ P2) Xk(q, -q + Pi + P2, Ps, Pa) 
ti2zin-mh exp(-(g - Pl - p 2 f/k 2 ) (exp(-g 2 /A 2 ) - exp(-g 2 /A; 2 )) 

+ (q^Pi+P2~q) (5.9) 



12 For this choice Rk(q) does not uniformly approach zero for k — > 0. The divergence for Rk(q 
0) is of no relevance in the present context 
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In the model under consideration this flow equation for A& will be integrated with 
a boundary condition for A& at k = A ~ 1.5 GeV, which is assumed to be generated 
after the gluons in QCD have been integrated out. Its form is motivated by the sum 
of one gluon exchange and a linearly rising potential, and it is assumed to depend 
on the variable t = (pi — p^) 2 only. After Fierz-transforming the one-gluon exchange 
diagram in spin and colour space and extracting the contribution proportional to 
the spin, colour, and flavour structure (|5.3|), this boundary condition reads 



Aa(Pi,P2,P3,P4) = + ~jT + D (t) ( 5 - 10 ) 

Here a s is the strong gauge-coupling constant which we take a s ~ .3, for the string 
tension we use A ~ .18 GeV 2 , and D{t) is a distribution with support at t — only, 
which is related to a constant in the potential [27] and will be given implicitly later. 

Solving the evolution equation fl5.9|) numerically - for details see the next section 
- we see indeed a pole-like structure as in eq. (1.5) appearing in A&. At the scale k v , 
where the collective fields are introduced, two conditions should be fulfilled: first, 
At should approximately factorize in the form 

A^(pi,P2,P3,P4) = g(PuP2)G(s)g(p 3} p 4 ) (5.11) 

Second, A^^ and hence T±u should not yet have become extremely large: T 4 appears 
on the right hand sides of the evolution equations for the higher N point functions 
r^v with N > 4. If T 4 becomes large, these higher N point functions necessarily 
become large as well, and their neglect becomes a very questionable approximation 
in this regime. (Generally, the introduction of the here neglected suppression factor 
fk of eq. (4.5) and the use of eq. (4.6) at some scale Ak v > avoid a possible 
conflict between these two conditions.) Within the present model we find, however, 
that both conditions can be met simultaneously to a reasonable extend at a scale 
k = k^ ~ .63 GeV (see sect. (6) and fig. (4)). 

As a result of eq. (5.11) can be written as 



4,Av 




(27r) 5 (pi + p 2 - p 3 - Pa) 



{[^(-Pl)^(P2)]^(Pl,P2)G( S )^(-p 4 ,-p 3 )[^(P4)^ a (-P3)] 
-[^a(-Pl)75^(P2)]^(pi,P 2 )G'(s)g(-p4,-p3)[^(p4)75^(-P3)]} 

+r 4 , v (5.12) 

Now we perform the step to switch to the effective average action for quarks and 
mesons Tk v fi[ l l J , Since we only consider the quark contributions to the evolution 

of the scalar part of the effective average action, there is no advantage to introduce a 
nonvanishing infrared cutoff k for the collective fields a in the present approximation. 
We choose composite operators which read in momentum space 

/d^p — ■ 

Oi 5)6 [V>,<M = - / 4^9(p,q-p)^i(-ph^t(q-p) (5-13) 
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where the real function g(pi,p2) is normalized 

2(0,0) = 1 (5.14) 
such that G(s) has dimensions ~ 1/s. Inserting eq. ( |5.13| ) in (4.6), one obtains 



+\ J ^{O a fe (g)G(g 2 )0 6 «(-g) +0^(g)G(g 2 )0 ( t(-g) 

-2~a\{q)O a \-q) - 2^\{q)0^ \-q)}. (5.15) 
It is convenient to define a by 

a = -(<r + J), 

a (5) = Ua-J) (5.16) 

such that the field a a b transforms under chiral flavour transformations Ul(N) x 
Uji(N) as an (N, N) representation 

5a% = i(6 a v c - 9% c )a\ ~ io* e W b + 6 A c b ) (5.17) 



After using T ktp [ifj , ip] = T 2 ,k tp [ip,ip] + T Atkip [ifj,tp] and eq. (|5.12Q for r 4)Av [V> ip], the 



terms quartic in ip in Ji(),ip, <j) cancel and we are left with 

rU^M = / -0rM^^ + l^r^-'iiVaiq) (5.18) 

-0^9(~q,Q -pWa(q)[\(l + 7fcK(p) + \(1 ~ T^ttC-jMfa ~ p)} 

We have neglected here the term F' Ak in eq. ( |5.12| ) which will be dropped in the 
following. 

Except for the momentum dependence of g the effective average action for com- 
posite fields IY describes a standard fermion-scalar theory with Yukawa coupling 
but without scalar self-interactions. 

In order to compute r ,o[0, a], we use the flow equation for I\ [0, a] in the form 
of the fermionic version ( |C.5| ) of eq. (2.23) 

dtT kfi [0,*] = - Tr {^( r St =0 + ^)" 1 } (5-19) 

(2) 

We observe that T k L^ =0 does not mix the fermionic and bosonic parts and we 
therefore only need to consider the second functional derivative with respect to 
the fermions. Apart from the inverse free quark propagator the only contribution 
to r fc U=o originates from the Yukawa coupling ~ aipifj, which we assume to be 
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independent of k in the form specified by ( [5.181 ). In this approximation the relevant 

,(2) 

fc,o|y>=o 



(2) 

part of oii/)=o becomes independent of k and eq. (|5.19|) is easily integrated. The 



result is 

r*,o[0, a] = -Tr In ( ffj^ ) + ^CT 1 ^ (5.20) 



where we implemented the boundary condition from fl5.18|) 



r fcv ,o(0,a) = ^G-V (5.21) 

in agreement with eq. (4.6). 

The effective potential for the field a or its effective wave function normalization 
can now easily be obtained from eq. ( |5.20| ) for k = 0: In the case of the effective po- 
tential I^ 2 ) on the r.h.s. of eq. ( |5.20| ) has to be evaluated for constant configurations 



of the field a, and the momentum dependent collective field propagator G~ l (q 2 ) is 
only needed at q 2 = 0. In the case of the effective wave function normalization Zk 
contributions arise both from an expansion of the Trln term up to second order in 
the external momentum (0(q 2 )), and from a corresponding expansion of G~ 1 (q 2 ). 
Note that this second contribution is already present at the scale k v according to 
eq. ( |5.21D ; thus an often used " compositeness condition" Zk = [fM], |25| does not 



hold at k = k v . 

As we will see in the next section, this effective potential exhibits indeed a non- 
trivial minimum for a real diagonal a. This expectation value is invariant under the 
vectorlike flavour symmetry, but carries nonvanishing axial charges. Our model 
therefore describes dynamical chiral symmetry breaking in this QCD-motivated 
fermionic theory. 

6 Numerical computations and results 

In this section we describe in more detail the procedure outlined in the last section, 
and we give an account of the results of the numerical solution of the evolution 
equation. The technique for integrating numerically the flow equation ( |5.9| ) follows 
the one developed in 0. The function X k can only depend on six independent 
Lorentz-invariant products of the momenta, which we denote by 

s = (P2 + pi) 2 = (p 3 + Pa) 2 , t=(pi- p-if = (P2 - Pif 

Vi = p 2 1: v 2 = p 2 , v 3 = p 2 , f 4 = p\. (6.1) 

Then we switch from the function Xk{pi) to its Laplace transform with respect to 
the Lorentz five invariants t, Vi,V2,V3,V4'. 

POO 

X k (Pi)= / dhdhdlidlACkisMe-*- 1 ™- 1 ™- 1 ™- 1 ** (6.2) 
Jo 

After inserting Q6.2Q into (|5.9|) and using 



2 2 

"X2- _ e -§r r i/k 2 



q 2 Ji/A 2 



dae~ aq2 (6.3) 
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the flow equation for Ck(s,k) can easily be derived, since the d A q integration in 
( |5.9| ) becomes Gaussian. The dependence of Ck(s,k) on the variables U is then 
discretized, thus Ck(s,li) becomes a function living on a five- dimensional lattice. 
The integration of the flow equation corresponds to an algorithm for updating the 
lattice for each step from k 2 to k 2 + Ak 2 . The boundary condition for C\(s, /j), with 
A = 1.5 GeV, corresponding to (|5.10|) reads 



C A (s, h) = C{l t )8{h)5{l 2 )5{h)5{h) (6.4) 

with 

C(l t ) = 8n\l t + 2ixa s (6.5) 

As discussed in 0, a vanishing of C{l t ) for l t > h max corresponds to a vanishing 
of the potential V(r) in ordinary space for r ^> (/t ma J 1//2 , and is related to the 
distribution D{t) in eq. Q5.10|) . With the help of the formulas in a direct relation 
between l tmax and the depth d of the linear part Ar of the potential can be derived: 

l t max = ^ (6.6) 

This finite value for h max specifies implicitly the distribution D(t) in (5.10). We will 
use (6.6) for l tmax with d = 1 GeV. 

With these methods and parameters, we compute Ck{s,li) on a lattice with 
50 points in the l t direction, and 10 points in the h,l2,h and I4 directions. For 
k = k v ~ .63 GeV we observe that Ck v (s, U) approaches a factorized form 

C kv (s,k)~6(l t )g(h,l 2 )G(s)g(h,h) (6.7) 

and G(s) becomes large for s — > 0. In fig. 4 we show, as a function of the scale 
k,G(s = 0) as a full line. If C(s,li) can be written as a product as on the r.h.s. 
of (6.7), a measure / defined by / = [C(s, l\ = 1, Z 2 = • C{s,l\ = OJ2 = 
0, I [C(s, h = l,l 2 = 0, lj) ■ C(s, h = 0, l 2 = 1, should satisfy / ~ 1. In fig. 4 
we also show this measure / as a dotted line. Thus for k — > k^ \ k indeed develops 
the pole-like characteristic for a bound state and we can perform the transition to 
a description in terms of a linear cx-model coupled to quarks. 

The next step involves the computation of the effective scalar potential from 
( |5.20| ). In the present model contains the free fermionic propagator and the 
Yukawa coupling to a a b as in eq. ( |5.18| ): 

(T( 2 \ q , q ')r b = i5 a b (2^5\q-q') (6.8) 

+^(-9, q'){A(Q - 90 + (^) a b ^ - 9) + «(Q - 9') - A(?' - 9))7 5 } 

In addition is diagonal in colour space. We will be interested in the effective 
potential for the field a b a . This obtains for a spatially homogeneous configuration 
which we will take diagonal and real: 

* h a(p) = ^a(2vr)V(p) (6.9) 
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With Q51j| ) for r^ 2 ) and the form (|5.6|) of the equation ( |5.20| ) for the effective 
potential Vfc(cr) becomes 



with 



(a k) q 2 

Pp ' ^ = (exp(-g 2 /A 2 ) - exp(-g 2 /fc 2 )) 2 (6,11) 

We observe, in particular, the fc-dependence of the mass term at the origin; 
neglecting the g-dependence of g we have 

^ - ^ - ^« - ■ ( (x^£^i) - < 6 - i2 > 

which becomes negative for small enough k and G _1 (0). This induces spontaneous 
breaking of the chiral symmetry. For the model under consideration and the choice 
of parameters discussed before, the resulting effective potential exhibits indeed a 
nontrivial minimum (independent of Nf). In fig. 5 we show 14 (<r) for k = .63, .27, 
.18 and GeV. Numerically the minimum of Vo(cx) is found to be at 

ct = -18 GeV (6.13) 



Due to our convention ( |5.14| ) for g(0, 0), and the aipip coupling present in T[ip, ip, a], 



this number can directly be interpreted as a constituent quark mass. 

In order to compute the chiral condensate < tptp > from the vacuum expectation 
value o"o we first Fourier-transform the operator (|5.12 ) 



O a \x) =-iJ d'yd A z^{x + y)${x + z)g(y, z) 

9{y,z) = J ^0^exvi(qy-q^)g(- q , q ') (6.14) 

For constant g(—q,q') = 1 the operator O a b reduces to ip l a (x)ip/ \x) (up to the 
phase factor which appears as a consequence of our Euclidean signature convention). 
The smoothening due to the wave function g(y, z) gives a well-defined regularized 
meaning^ to the chiral condensate < ipip > as the expectation value of O 

< ^ >= _L < Tr 0(x) > (6.15) 
Nf 

The exact relation of our field-theoretical definition of < tfjip > to more phenomeno- 
logical definitions in the context of QCD sum rules remains to be established. We 
emphasize that our formalism offers, in principle, the possibility to provide well- 
defined field-theoretical definitions for the concepts underlying QCD sum rules or 

13 An additional smoothening will be related to the appearance of fk in the definition of O 
(4.4) in a more accurate treatment. 
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vacuum correlators. The relation between the vev of the collective field <J a b (x) and 
the expectation value of the composite operator O/ is given by eq. (3.11) 

< a a h >= (5(0) < O a b > . (6.16) 

Therefore the relation between o"o and the QCD quark condensate reads 

a = G(0) < > . (6.17) 

With (5(0) = 33.2 GeV" 2 we thus find 

< Tpip >= (175 MeV) 3 . (6.18) 

It should be clear that the collective fields a b a contain the massless pionic degrees of 
freedom, once they are expanded around the vev fl6.9|). We can also obtain the pion 
decay constant f n , if we expand the effective action r OiO (0, a) up to second order in 
derivatives acting on a. Let us denote the corresponding kinetic terms in I\o by 

With this normalization we have 

U = \J 2Zk=o ■ o"o (6.20) 

As stated below eq. (|5.21|) , Z gets contributions from both terms on the r.h.s. of eq. 
( p.2U| ). In the model under consideration we actually find a dominant contribution 
.825 arising from the term involving together with a contribution .017 from 

the fermionic loop we thus obtain 

f w = 234 MeV (6.21) 

which has to be compared with the experimental value f n = 93 MeV. We observe 
a substantial discrepancy, but expect that the dominant contribution arising from 
the term involving G~ l is particularly sensitive to the rather crude approximations 
employed in this paper. 

7 Conclusions 

We have presented in this paper a formalism to describe bound states and con- 
densates of composite operators in terms of an effective average action. The scale 
dependence of this effective action is determined by an exact nonperturbative evo- 
lution equation. The properties of the vacuum and excitations around it can be 
infered from the solution of the flow equation for k — > 0. Such a solution allows to 
extrapolate from the short distance physics to the long distance physics. We have 
deviced a detailed prescription how a change of relevant variables can be performed 
in the course of this evolution. 
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The practical use of the exact flow equation is closely related to the existence 
of a suitable truncation scheme which permits to solve it approximately. We have 
demonstrated the existence of nonperturbative approximation schemes for a QCD- 
motivated quark model. Our computations in sect. 5 and 6 may be viewed as 
first steps in a systematic expansion in 1PI Green functions with a fixed number of 
external legs for the quark fields. Here the truncation neglects many-quark-operators 
but retains the full momentum dependence of the lower 1PI vertices. On the other 
hand the truncation for the meson part of the effective action may be developed 
as an expansion in momenta (or appropriate functions of momenta). This retains 
the full field dependence on the zero momentum modes and can therefore give a 
detailed picture of the meson potential and kinetic term. Although in the present 
paper we use very rough approximations we have demonstrated that our formalism 
provides not only a formal tool to describe phenomena as complex as dynamical 
chiral symmetry breaking, but that suitable truncations of the fc-dependent effective 
action allow also to compute phenomenologically interesting numbers. 

Here, however, no effort was made to discuss systematically the dependence of 
the results on various parameters entering the approximations, like A, k v or the 
string tension A. We defer such an investigation to future work, where less drastic 
approximations to the effective action will be made. Such extensions do actually 
not necessarily involve an enormous amount of additional computational effort. The 
following parts of such a programme are certainly feasable and have partially already 
been performed: In the case of the four point function more spin, flavour and 
colour structures can be taken into account. After the introduction of the collective 
fields the fc-dependence of the Yukawa-like couplings and the quark wave function 
normalisation do not have to be neglected. 

The meson self-interactions can be included on the r.h.s. of the flow equation, 
using the picture with infrared cutoff for both fundamental and composite fields de- 
veloped in sect. 4. The transition from the description in terms of only fundamental 
fields to a language with collective fields can be smoothened by an appropriate sep- 
aration of the bound-state scale from the scales k v , Ak v as discussed in sect. 5. 
Quark mass terms can be included without any conceptual difficulties. We are also 
free to include the p-mesons as additional collective fields. Finally, we may replace 
the "initial conditions" for the four quark vertex at the ultraviolet cutoff scale A 
by a solution of the evolution equation for QCD, i.e. the coupled system of quarks 
and gluons. The practical applicability of the effective average action for nonper- 



turbative problems in gauge theories has already been demonstrated [26]. Although 



this program still needs many steps it seems not impossible to us that it may finally 
lead to a reliable computation of chiral condensates and the meson spectrum and 
interactions. 
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Appendix A: Green functions for collective fields 
and composite operators 

The generating functional T [<p, a] contains all information about Green functions for 
the collective field a, and hence as well as for the composite operator O. The relation 
between the two sets of Green functions is given in this section. The generating 
functional for the 1PI Green functions for the collective field obtains from r [<^,cr] 
as 

T[a] = T [ip [a],a] 

5( P \<po 

In particular, the minimum of Tq 



6To 

5a \cp,a 
5T 





(A.2) 



5< P \w 

determines the expectation values of x an d 0\x] 

<x a > = ^ 

<d[xY> = a*. (A.3) 

The higher correlations for the composite operator 0[x] are not directly given by the 
functional derivatives of T[a] with respect to a but they are simply related to them. 
From fl2.1|) it follows immediately that the generating functional for the connected 



Green functions for O is 

W[J, K\ = W [J, K\ - ]^K ] GK. (A.4) 

Comparing the Legendre transforms of W and Wo one finds that the generating 
functional for the 1PI Green functions for O is given by T 

f = -W + JV + K ] a 

= -W + ^K j GK + Jfy + J&o - K ] GK 

= To-^GK (A.5) 



Here we have used 



dW 
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The generating functional for tp and 0[<p] is therefore obtained from To be straight- 
forward algebraic manipulations 



t[<p,a] 



2 5a* 



Ha] 



a 1 



° ~^5al 



(A.7) 



(A. 



As expected, the difference between F and T is irrelevant for the expectation values 
but it accounts for ( 2.14 ) and generalizations to higher 1PI vertices. The 



Green functions for O obtain from the Green functions for O by simple rescaling 



Appendix B: Exact evolution equation restricted 
to composite fields 



If (fo depends on a one has to correct fl4.3|) according to 

da = + W 1 ^ 1 ^ = ° (R1) 

and use 

d 2 T[a] d 2 T r . . . d 2 T . . . .dtp ,„ . 

since d/da = d/da + (dp /da)d/dp. The evolution equation ( |3.7| ) expressing 
dTj,[a\/dt in terms of r| [<r] remains nevertheless exact. This can be seen most 
easily by using an implicit functional integral representation for T~ k [a] 

exp -V~ k [a] 

= J VxVpexp - is[ X ] + \g~\p - GO[ X }) 2 - ^l(p - a) + A k S[p - a]} 

= jvp'e^-\s efi \p + p'\- 5 ^p' + K- k S\p']\. (B.3) 

We note the complete analogy with the functional integral representation [16 of the 
effective average action for the fundamental fields 

exp-r*M= IPx'exp-j^ + xl-^^x' + A^x']}. (B.4) 

The exact evolution equation can be infered from this analogy by observing that the 
precise functional form of S e ff needs not to be known. Alternatively, one can derive 
the evolution equation starting directly from ( |B.4| ) and following the procedure of 



ref. 111. One arrives at 



(2) 

with r~ now denoting second functional derivatives with respect to a. 
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Appendix C: Effective average action with fermions 



In this short section we present the generalization of our formalism for fermions. We 
write the fermionic infrared cutoff term as 

A k S F = Q fj a (R kF ) a ( flf. (CI) 

Here rj are Grassmann variables which are either Dirac spinors or Weyl spinors 
in even dimensions, or Majorana or Major ana- Weyl spinors in those dimensions 
where this is possible . (We use the conventions for Euclidean spinors of ref. f2~2fl .) 
The factor | appears only if r] and fj are related, as for the case of Majorana and 
Majorana- Weyl spinors or a representation of Weyl spinors in the form of 2^ com- 
ponent spinors with identified components. For Dirac spinors or the standard 2a _1 
component Weyl spinors this factor is absent. The formalism of ref. j| applies 
identically for fermions except for one important minus sign. First we observe that 
the infrared cutoff does not mix bosons and fermions. The evolution of the effec- 
tive average action dT^/dt obtains therefore separate contributions from bosons and 
fermions. Both for the bosonic and fermionic contribution one arrives at the identity 

|r fc = \j t {R kB )%{< xlx p > - < xl >< >) 

+ Q ^f)V< W >-<fh><rf >)• (C2) 

The minus sign for the fermionic contribution appears due to the anticommutation 
properties of the Grassmann variables 

< f} a rf > - < fj a >< if >= -(< vff\ a >-</>< f] a >) = -((r k 2) + R k ) FF f a 

(C.3) 

(2) 

Here T k and R k are considered as enlarged matrices acting both on bosonic and 
fermionic degrees of freedom, i.e. 



f F (2)\/3 
\ L k JFFa 



/ F (2K/3 
\ L k )BFo 



\ L k I FBa 



5 2 r, 6 2 T k 



5 2 T k 



S 2 T k 



(If } )V = -f^~ a (C.4) 

(We observe that the second derivatives ( |C.4| ) are sufficient to describe T k only if 
X, X* an d ViV are n °t independent, e.g. for real scalar fields and Majorana spinors. 
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For Dirac spinors r\ and fj lead to independent entries in the matrix of second func- 
tional derivatives and one needs, for example, 5 2 Tk/5x%Sfj a .) The inclusion of possi- 
ble fermionic composite operators is straightforward and our final evolution equation 
is 




(C.5) 



The first trace effectively runs only over bosonic degrees of freedom and the second 
over fermionic ones. The second functional derivative 1^ involves both fundamental 
and composite fields. The evolution with respect to k is similar (cf. (|3.7|) ). 



Figure Captions 

Fig. 1: Vertex joining the collective field p to two fundamental fields x as present 
in the action Sk[x, p] of eq. (3.3). 

Fig. 2: Form of the term ~ 0^[%]GO[x] in the action Sk[x, p] of eq. (3.3). 

Fig. 3: Diagrammatic representation of the evolution equation (5.9) for the four 
point function X k . 

Fig. 4: Full line: The collective field propagator G(s), at s = 0, as a function of the 
scale k, in units of 100 GeV~ 2 . This result has been obtained from a numerical 
computation of C& as described in sect. 6, and using the relation (6.7) between 
Cfc and G. 

Dotted line: The measure /, as defined below eq. (6.7), as a function of the 
scale k. The approach of / towards 1 indicates that the momentum dependence 
of the four point function Xk or its Laplace transform Ck indeed factorizes as 
used in eq. (6.7). 

The scale k v , where the collective field is introduced, is chosen at .63 GeV. 

Fig. 5: Effective potential for the collective field a, as given by eq.(6.10), for 
different values of the infrared cutoff k. The dotted lines correspond to, from 
inside to outside, k = .63 GeV, k = .27 GeV and k = .18 GeV. The full line 
is the final result for k = 0. 
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